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Abstract Pulsar braking indices offer insight into the 
physics that underlies pulsar spin-down. Only five brak- 
ing indices have been measured via phase-coherent tim- 
ing; all measured values are less than 3, the value ex- 
pected fi'om magnetic dipole radiation. Here we present 
new measurements for three of the five pulsar braking 
indices, obtained with phase-coherent timing for 
PSRs J1846-0258 {n = 2.65 ± 0.01), B1509-58 (n = 
2.839±0.001) and B0540-69 {n = 2.140±0.009). We dis- 
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cuss the implications of these results and possible phys- 
ical explanations for them. 
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1 Introduction 

A very commonly assumed model for pulsar spin-down 
posits that 



(1) 



where v is the pulse frequency, i> is the frequency deriva- 
tive, is a constant and n is the braking index. The 
braking index is then given by 



1/1/ 



(2) 



The braking index provides insight into the physics that 
drives pulsar spin-down. Typically, it is assumed that 
magnetic dipole radi ation underlies pulsar evolut ion, re- 
sulting in n = 3 fe.g. [Manchester fc Tavloij|l977l ). How- 
ever, other processes could, in principle, cause the pul- 
sar to radiate and would result in different values for 
n and K. For example, a pulsar spun down entirely 
by the loss of relativis tic particles would have n = 1 
(jMichel fc TuckeI^ll969^ . A pulsar losing energy via grav- 
itational radiation or q uadrupole magnetic radiatio n would 
spin down with n = 5 ( Blandford fc Romanilll988[ ). 

Braking indices have proven difficult to measure. To 
date, only six have been reported even though more than 
1600 pulsars are known. Evidently, the pulsar properties 
necessary for a measurement of n are rare; the pulsar 
must: spin down quickly; experience few, small, and rel- 
atively infrequent glitches; and be relatively uncontam- 
inated by timing noise, a low-frequency stochastic pro- 
cess superposed on the deterministic spin-down of the 
pulsar. The youngest pulsars, of which the Crab pulsar 
is the most famous example, uniquely possess these three 
qualities. The reason then for the paucity of measured 
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values of n is a direct consequence of the relative rarity 
of very young pulsars, i.e. those with characteristic ages 
on the order of 1 kyr, where characteristic age is defined 

as 



P 
2P 



2v 



(3) 



Of the six pulsars with measured n, five were ob- 
tained via phase-coherent timing. All five of these pul- 
sars: PSRs J1846-0258, B0531-I-21 (the Crab pulsar, 
n = 2.51 ± 0.01), B1509-58, J1119-6127 (n = 2.91 ± 
0.05), and B0540— 69, have characteristic ages less than 



3 Observations 

In this paper we describe observations of three young 
pulsars PSRs J1846-0258, B1509-58, and B0540-69 
taken with the Proportional Counter Array (PCA) on 
board the Rossi X-ray Timing Explorer (RXTE). The 
PCA consists of five coUimated xenon/methane multian- 
ode proportional counter units (PCUs). The PCA oper- 
ates in the 2-60 keV energy range, has an effective area of 
~6500cm^ and has a 1 degree field of view. While RXTE 
has no imaging capability, it has e xcellent time resolu- 
tion of - 1 /iS (|Jahoda et al.lll996l ). This makes RXTE 



2 kyr (Livingstone et al.l2006l:[Lvne et al.lll993l:lLiyingstoneideallf or observing young, rapidly rotating pulsars. 
l2005bl : ICamilo et al.ll2000l : 'Livings tone et al.ll2005al ). The 
sixth measurement, that of the Vela pulsar (n = 1.4 ± 
0.2), could not be o btained with pha se-coherent timing 
due to large glitches ( Lvne et al]ll996t ). Timing noise and 
large glitches begin to seriously contaminate measure- 
ments of n when pulsars have characteristic ag es ^ 5 kyr 
( McKenna fc Lvnd 1 1990l : [Marshall et al.ll2004ti . 

In this paper we report on long-term Rossi X-ray 
Timing Explorer observations of three young pulsars, 
PSRs J1846-0258, B1509-58 and B0540-69. We present 
braking index measurements for each of these pulsars ob- 
tained via phase-coherent timing. 



2 Phase-coherent pulsar timing 

The most accurate method of extracting pulsar timing 
parameters is phase-coherent timing, that is, account- 
ing for every turn of the pulsar. Pulse times of arrival 
(TOAs) are measured and fitted to a Taylor expansion 
of pulse phase, 4> at time t given by 



2 D 



where subscript denotes a parameter at the reference 
epoch, tg. TOAs and initial spin parameters are input to 
pulse timing software (e.g. TEMPC0) and refined spin 
parameters and timing residuals are output. 

The existence of timing noise and glitches in young 
pulsars is well known to contaminate the measurement of 
deterministic spin parameters. Though powerful, a fully 
phase-coherent timing solution can be sensitive to these 
contaminants. In such partially coherent method 

may be employed. In this case, local phase-coherent mea- 
surements of v, v and possibly i> are made. Though the 
efi'ects of timing noise cannot be eliminated, the noise 
component is more readily identified and separated from 
the deterministic component of the spin-down, as shown 
in Section [5l In addition, glitches can be easier to iden- 
tify with this method, as will be shown in Section [6] of 
this paper. 

^ http:/ /www. atiif.csiro.au/research/pulsar/tempo 



Observations of PSR B1509— 58 were taken in 
"GoodXenonWithPropane" mode, while observations of 
the other two sources were taken in "GoodXenon" mode. 
Both modes record the photon arrival time with l/xs- 
resolution and photon energy with 256-channel resolu- 
tion. The number of PCUs active during an observation 
varies, but is typically three. For PSRs B1509— 58 and 
J1846— 0258, which have relatively hard spectra, all three 
Xenon layers and photons with energies ranging from 2- 
60 keV were used, while for the softer spectrum source, 
PSR B0540— 69, only the top Xenon layer and photons 
with energies ranging from 2-18 keV were used. Further 
details of X-r ay and radio observ a tions o f PSR B1509— 58 
are given in [Livingston e et al" 1 (|2005bh and references 
therein. Details of RXTE observations of PSR B0540-69 

can be found in lLivingstone et aTl ( 2005a( ) whil e details of 

PSR J 1846— 0258 observations can be found in lLivingstone et al.l 
(I2006D and references therein. 

Data were reduced using standard FITS tools as well 
as in-house software developed for analyzing RXTE data 
for pulsar timing. Data from different PCUs were merged 
and binned at (1/1024) s resolution. Photon arrival times 
were corrected to barycentric dynamical time (TDB) at 
the solar system barycenter using the J2000 source po- 
sitions and the JPL DE200 solar system ephemeris. 

For PSRs B0540-69 and J1846-0258, initial ephemerides 
were found by performing periodograms on observations 
to determine values of v. Several values of v were fit- 
ted with a linear least squares fit to determine an initial 
value of V. These initial values were then used as input 
to a Taylor expansion of TOAs to determine more accu- 
rate parameters (Eq.|3]). PSR B1509— 58 has a previously 
determined ephemeris from radio timing data obtained 
with the Molonglo Observat ory Synthesis Tele scope and 
the Parkes Radio Telescope ( Kaspi et al.lll994l) . We were 
able to extend that fit with 7.6 yr of RXTE data, remov- 
ing a constant, but not well determined, offset between 
radio and X-ray TOAs. 



4 PSR J1846-0258 

PSR J1846-0258 is a very young pulsar (r^ = 723 yr) lo- 
cated at the center of the supernova remnant Kesteven 
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Table 1 Spin parameters for PSR J1846-0258 



Parameter 



First solution 



Second solution 



Dates (MJD) 
Epoch (MJD) 

V (s-i) 

V (10-"s-2) 

V (10-2is-3) 
Braking Index (n) 
Glitch epoch (MJD) 
Avjv 

Avjv 



51574.2 - 52837.4 
52064.0 

3.0782148166(9) 

-6.71563(1) 

3.87(2) 

2.64(1) 

52210(10) 

2.5(2) X lO"'' 

9.3(1) X 10"'' 



53578.6 



52915 
53404.0 
3.070458592(1) 
-6.67793(5) 
3.89(4) 
2.68(3) 



2000 



Epoch (Years) 
2001 2002 



2003 



200 



100 



I I I \ ^ ^ ^ i p 
RMS = 0.146 periods 



T" 



75. It has a relatively long spin period of 324 ms and a 
large magnetic fielcfl of B - 5 x lO^^ G. PSR J1846-0258 
has been ob served with RXTE fo r 6.3 yr since its discov- 
ery in 1999 (|Gotthelf et al.ll2000l) . 

Using our initial ephemeris we obtained a 
phase-coherent timing solution valid over a 3.5 yr inter- 
val in the range MJD 51574-52837. Three spin parame- 
ters (y^ V and v) were required by the fit. In addition, 
we discovered a small glitch near MJD 52210±10. The 
fitted glitch parameters are Avjv = 2.5(2) x 10~^ and 
Ai>/i> = 9.3(1) x 10"^. Note that these and aU other 
quoted uncertainties are 68% confidence intervals, unless 
otherwise indicated. The wide spacing of data near the 
glitch prevent the detection of any short-timescale glitch 
recovery. Timing residuals are shown in Figure[TJ The top 
panel of Figure [T] shows residuals with i>, i) and glitch 
parameters fitted. The residuals clearly show systemat- 
ics due to timing noise and possibly unmodeled glitch 
recovery. In order to minimize contamination of long- 
term timing parameters, we fitted additional frequency 
derivatives to render the residuals consistent with Gaus- 
sian distributed residuals (a pro cess known as 'wh iten- 
ing' residuals. See, for example, iKaspi et al.lll994l ). For 
this timing solution, a total of eight frequency derivatives 
were fitted, shown in the bottom panel of Figure [TJ The 
braking index resulting from this 'whitened' timing solu- 
tion is n = 2.64 ± 0.01. Complete spin-down parameters 
are given in Table [1] 

Phase was lost over a 78-day gap in the data near 
MJD 52837, indicated by the fact that a timing solution 
attempting to connect over this gap fails to predict the 
pulse frequency at previous epochs. The loss of phase is 
likely due to timing noise or a second glitch. However, 
these two possibilities cannot be distinguished due to the 
relatively long gap in the data set. 

A second phase-coherent solution was obtained for 
1.8 yr from MJD 52915-53579 with v, v and v. Timing 
residuals with these three parameters fit are shown in the 
top panel of Figure [21 Again, systematics due to timing 
noise and/or glitch recovery remain in these residuals. To 
'whiten' timing residuals, five total frequency derivatives 
were fitted from the data, shown in the bottom panel 
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Fig. 1 Phase-coherent X-ray timing analysis of the young 
pulsar PSR J1846— 0258 spanni ng a 3.5-yr interv al in t he 
range MJD 51574-52837 (after ILivingstone et all ([200^ )'). 
Top panel: Residuals with v, v as well as glitch parameters 
Av and Ail fitted. The glitch epoch, MJD 52210 is indicated 
by the arrow. Bottom panel: Residuals with glitch parame- 
ters and eight frequency derivatives in total fitted to render 
the residuals consistent with Gaussian noise. 



of Figure [2] Complete spin parameters for this timing 
solution are given in Table [TJ 

Quoting the average value of n from the two inde- 
pendent timing solutions, which are in agreement, gives 
n — 2.65 ± 0.01. As is the case for all measured values of 
n, this value is significantly less than 3, the value consis- 
tent with spin-down via magnetic dipole radiation. This 
implies that some other physical process must contribute 
to the spin-down of all of these pulsars. 

This measur ement of n for PSR J184 6— 0258 increases 
its age estimate ( Livingstone et aO2006[ ). The commonly 
known characteristic age (Eq. [3]) implicitly assumes that 
n — "i. K more physical estimate can be made once n is 
known. The age then can be estimated as 



1 



1 



Tc < 884 yr. 



(5) 



B EE 3.2 X 10^^(PP)^''^G 



The estimate is an upper limit since the initial spin 
frequency of the pulsar is not known. The upper limit 
approaches an equality when the pulsar is born spin- 
ning much faster than its present spin frequency. Given 
the long period of the pulsar and the estimated initial 
spin period distribution (e.g. iFaucher-Giguere fc Kaspil 
I2006D . the latter is likely to have occurred. This age es- 
timate for PSR J1846— 0258 is less than the known age 
of the Crab pulsar of 952 yr. 
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Fig. 2 Phase-coherent X-ray timing analysis of 
PSR J1846— 0258 spann ing an 1.8-yr i nterva l from 
MJD 52915-53579 (after ILivingstone et al.l IpOT )). Top 
panel: Residuals with f, !>, i> fitted. Bottom panel: Residuals 
with five frequency derivatives total fitted to render the 
residuals consistent with Gaussian noise. 
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Fig. 3 Timing residuals for PSR B1509-58. Radio TOAs 
are sh own eis dots and the X-ray TOAs are shown as crosses 
(after [Tivingstone et al.l (|2005a )). The top panel has pulse 
frequency and three frequency derivatives removed, the mid- 
dle panel also has the fourth frequency derivative removed, 
and the bottom panel shows residuals after the removal of 
five frequency derivatives. 



5 PSR B1509-58 



The young pulsar PSR B1509— 58 was discovered in 1982 
and has been observed regularly ever since, first with 
radio telescopes such as the Moloii^lo O bservatory Syn- 
thesis Telescope (Manchester et al.' 1985) and the Parkes 
Radio Observ atory (.Kaspi et al. .l994,). and more re cently 



with RXTE (Roti 



y iKaspi et a.L.L\j\jy. and m ore recently 
^ 120041 : ILivingstone et al.l l2005bt) . We 



phase-connected all 21.3 yr of available radio and X-ray 
timing data to determine the braking index. Timing resid- 
uals are shown in Figure [31 The top panel shows the 
timing residuals with v and three frequency derivatives 
fitted; the middle panel shows residuals with the fourth 
frequency derivative also fitted; the bottom panel shows 
timing residuals with five frequency derivatives fitted. 
Remarkably for such a young pulsar, no glitches were 
detected in this time period. Spin parameters from this 
phase-coherent analysis are 
1^ = 6.633598804(3)s-i, u = -6.75801754(4) xlO^ 
V = 1.95671(2) X 10-21 at epoch MJD 49034.5. These 
parameters imply a braking index of n = 2.84209(3), 
though timing noise that could not be completely re- 
moved by fitting additional frequency derivatives con- 
tributes to a systematic uncertainty that is not included 
in the formal uncertainty quoted here. To solve this prob- 
lem, we performed a partially phase-coherent analysis by 
making independent measurements of n. 



'11 c 



Due to the large value of z> for this pulsar, a signif- 
icant measurement of n can be made in approximately 
2yr, without noticeable contamination of the measured 
spin parameters from timing noise. Thus, having over 20 
years of data allows 10 independent measurements of n, 
which are shown in Figure 21 No secular variation of n 
over 21.3 yr is seen, however, there is significant deviation 
from the average value of n = 2.839 ± 0.003. This uncer- 
tainty was determined by a 'bootstrap' analysis which is 
a robust method of determining uncertainties when the 
formal uncertainties are thought to underestima te the 
true v alues, i.e. due to the presence of timing noise (|Efronl 
I1979D . Note that this value is in agreement with that ob- 
tained with the fully phase-coherent timing solution, as 
well as the previou sly reported value of n = 2.837±0.001 
(jKaspi et al.lll994l ). The reduced is 15 for 9 degrees 
of freedom. This variation is likely due to the same tim- 
ing noise process that can be observed in timing residu- 
als. Here, the variation is at the ~ 1.5% level. A similar 
analysis has been performed for PSR J1846— 0258 where 
variations are seen t o be on the order of ^ 5%, though 
only at the 2a level ([Livi ngstone et ani2006l ) and for the 
Crab pulsar wher e variations are on the order of 0.5% 
(|L^ne!eEal][i99l). 
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46000 48000 50000 52000 
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Fig. 4 Braking index calcula ted at 10 epochs of '~2yr in 
length (after Livingsto ne et al.l ([2005b)). There is no statisti- 
cally significant secular change of 21.3 yr of data. The average 
value is 2.839 ± 0.003, in agree ment with the previously re- 
ported value (IKaspi et al.lll99^ and the value obtained from 
a phase- coherent analysis. The reduced x'^ value is 15 for 9 
degrees of freedom, suggesting contamination by timing noise. 



the 300 days of data available after t he glitch to mea- 
sure a braking index of n = 1.81 ±0.07. ICusumano et"al] 
(I2003D extended the data set to 4.6 yr and reported that 
no glitch occurred. In contrast to the previous value, they 
measured n = 2.125 ± 0.001. 

We re-examined all previously reported RXTE data 
and extended the data set by Syr in order to resolve 
the discrepant timing solutions and measure the true 
braking index for this source. We phase connected a 
total of 7.6 yr of data and found a small glitch near 
MJD 51335 with parameters A1//1/ ~ 1.4 x 10"^ and 
A vjij ^ 1.33 X 10~^, in agreement with those reported 
by IZhang et al. 1 (|200ll ). This glitch is very small, and is 
most easily seen by the change in v at the glitch epoch. 
Figure [5] shows 22 measurements of v obtained from in- 
dividual phase-coherent analyses, with the fitted glitch 
epoch indicated by an arrow. The slope of the line, that 
is, the second frequency derivative does not change sig- 
nificantly after the glitch (before v — 3.81(3) x 10~^^s~^, 
after v = 3.81(1) x lO^^^s""^). Uncertainties on i) were 



6 PSR B0540-69 

Located in the Large Magellenic Cloud, PSR B0540-69 
is commonly known as the 'Crab Twin', due to its simi- 
lar spin and nebular properties. For instance, its period 
of 50 ms and magnetic field _B 5 x 10^^ G are nearer 
to those of the Crab pulsar (P 33 ms, B ~ 4 x lO^^ q) 
than for any other pulsar. Due to its large distance, 
PSR B054 0-69 is very difficult t o detect in the radio 
waveband ( Manchester et al.1ll993f ). hence regular radio 
timing of this source is not practical. 

The lack of regular, long-term timing observations for 
this pulsar has led to conflicting values of n in the lit- 
erature; report ed values raiige from n = 1.81 ± 0.07 to 
n = 2 .74±0.01 (jZhang et al.ll200ll : l Ogelman fc Hasingeil 
I1990I ). Widely spaced timing observations greatly increases 
the risk of losing phase if a phase-coherent solution is at- 
tempted. If instead of a phase-coherent timing solution, 
measurements of frequency are obtained over widely spaced 
intervals, small glitches can easily be missed and the ef- 
fects of timing noise are difficult to discern. Two con- 
flicting values of n are of particular interest since they 
are based on ov e rlappi ng data from RXTE. 

[Zh ang et all (j200ll ) reported on 1.2 years of regular 
timing observations and found a small magnitude glitch 
at MJD 51325±45 with parameters Ai^/f = (1.90 ± 
0.04) X 10~9 and Ai>/i> = (8.5 ± 0.5) x 10"^. They used 



determined by a bootstrap analysis ( Efron|[l979l ) . We use 
the average to determine the braking index, found to be 
n = 2.140 ±0.009. 

In agreement with IZhang et al. I (|200lh . we report a 
small glitch near MJD 51335, though our value of n is 
signiflcantly larger. By phase-connecting only the same 
300 day subset of data that they used to measure n, we 
find n = 1.82 ± 0.01, in agreement with their result. The 
low value of n in this case appears to be the result of 
timing noise and/or glitch recovery contaminating the 
relatively short time baseline used to measure n. 

Our measured va l ue of n is 1.7a from that reported 
bv ICusumano et all (|2003l . n = 2.125 ± 0.001) though 
they do not report a glitch and their uncertainty does 
not account for the effects of timing noise. The rea- 
son for the agreements between our measured values 
is that their value of n was determined by two phase- 
coherent fits to the data, b e fore a nd after the glitch epoch 
reported by Zhaiig et al.l ( 200ll ). despite the fact that 
ICusumano et al. (loojreport no glitch. 



7 Implications and Physical Explanations for 
n < 3 

All measured values of n are less than 3, the value consis- 
tent with spin-down due solely to magnetic dipole radia- 
tion. This implies that an additional torque is contribut- 
ing to the spin-down of young pulsars. Also intriguing is 
the relatively wide range of measured values of n, shown 
in Table [2l A measurement of n immediately provides 
a correction to the age estimate of the pulsar given by 
the characteristic age, as shown with Equation [5l Com- 
parisons of the characteristic age and age estimated with 
n are also given in Table [2l Although the age estimate 
is always increased with a measurement of n < 3, it 
should be noted that the age estimates given are upper 
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Fig. 5 Measure ments of i>; the slope is i> (after 
ILivingstone et all (|M)5a)). The ghtch occurring near 
MJD 51342 is shown with an arrow. The pre-ghtch slope 
is i> = 3.81(3) X 10~^^s~'', while the post-glitch slope is 
u = 3.81(1) X lO'^^^s"^. The average of pre- and post-glitch n 
is 2.140 ±0.009. Measurement uncertainties are smaller than 
the points. 



Table 2 Braking index measurements via phase-coherent 
timing. Also given are the characteristic age, Tc and age esti- 
mate using n, r. 



Pulsar 


n 




T 


J1846-0258 


2.65(1) 


723 


884 


B0531+21 


2.51(1) 


1240 


1640 


B1509-58 


2.839(3) 


1550 


1690 


J1119-6127 


2.91(5) 


1610 


1680 


B0540-69 


2.140(9) 


1670 


2940 



limits since the initial spin frequencies are not known. 
The calculation of magnetic fields are also affected by a 
measurement of n < 3, since these are obtained assum- 
ing pure magnetic dipole radiation. Unfortunately, there 
is no simple formula to estimate the correction to the 
dipole magnetic field as there is for the age. Specific de- 
tails of the spin-down torque are required to uncover the 
true magnetic field of pulsars. 

There are several theories that attempt to explain the 
measurements of n < 3. One explanation is that the pul- 
sar's magnetic field grows or counter-aligns with the spin 
axis. This is equivalent to allowing the 'constant', K, in 
the si mple model of pulsar spin do wn (Eq.[T]) to vary with 
time ( Blandford &: RomanilllQ SSh. An advantage of this 



To date, the third frequenc y derivative has been 
measured only fo r the Crab pulsar (iLvne et al.lll993h and 
PSR B1509-58 (jLivingstone et all l2005tO hough nei- 
ther is known with sufficient precision to rule out the 
null hypothesis of a constant value of K. Timing noise 
and in the case of the Crab pulsar, glitches, may prevent 
a sufficiently precise measurement from ever being made. 

Another suggestion is that a fall-back disk formed 
from supernova material modulates the spin-down of young 
pulsars, providing a propeller torque in addition to the 
torque from magnetic dipole radiation. This would cause 
the pulsar to lose ener gy more quickly le ading to a mea- 
sured value 2 < 71 < 3 (jAlpar et al]l200lh . A difficulty in 
this model is that the disk must not suppress the pulsed 
radi o emission during the propeller phase (Mcnou et af] 
[200 ih . 

In recent years, much work has been done on mod- 
elling the pulsar magnetosphcre. Fully physical, three 
dimensional, time-dependent models of the pulsar mag- 
netosphcre are still some time away, however, significant 
progress has been made and there is some suggestion that 
n < 3 may be a natural re sult o f a plasma filled mag- 
netos phere (see, for example SDi1:kovskv..2005. : .TimokhinI 
120061 : IContopoulos k Spitkovskvll20oi L 



model is that it can be tested if precision me asurements 
of the third frequency derivative can be made (iBlandfordI 



The idea that plasma in the magnetosphcre affects 
the torque acting on a pulsar is gaining acceptance with 
the first observa tional evidence for th is having recently 
being presented. iKramer et al. 1 (I2006D show that 
PSR B1931-I-24, which has curious quasi-periodic nulling 
behaviour, spins down at different rates when it is 'on' or 
'off'. Specifically, the pulsar has a faster rate of spin down 
when it is observed in the radio waveband than when 
it goes undetected. This implies a connection between 
the radio emission mechanism and the spin-down torque. 
The interpretation presented by the authors is that the 
radio emission mechanism is only active when sufficient 
plasma is present in the magnetosphcre, and that this 
plasma exerts a torque on the pulsar, spinning it down 
faster than in the absence of plasma. If this is indeed the 
case, then all observable pulsars should have n < 3. 

iMelatosI ([1993) suggested that the solution to the 
n < 3 problem is related to the angle between the spin 
and magnetic axes, a, as well as to currents in the mag- 
netosphcre. Melatos postulates that the magnetosphcre 
can be considered to be split into two sections, an in- 
ner and outer magnetosphcre. The division occurs at 
the 'vacuum' radius, the location where particles are no 
longer confined to field lines. The inner magnetosphcre 
will then corotate with the neutron star and can be con- 
sidered part of the radius of the rotating dipole. However, 
since this radius is less than, but comparable in size to, 
that of the light cylinder, the dipole can no longer be 
treated as a point, but has some finite size. As a result, 
2 < n < 3 and n approaches 3 as a pulsar ages. This 
model is especially attractive because it provides an ex- 
planation for the large scatter in observed values of n, 
and provides a prediction for n given measured values of 
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z/, V and a. Given the large uncertainties on known val- 
ues of a, the model roughly agrees with measurements 
of n for PSRs B1509-58, B0540-69 and the Crab pul- 
sar. PSR J1119— 6127 does not appear to have a well 
determinable a. However, following the Melatos model, 
the measured valu e of n = 2.91 ± 0.05 predicts a range 
of 10° < a < 32° (jCrawford fc Ke im 2003). Our recent 
measurement of n for PSR J1846— 0258 allows a predic- 
tion of a = 8.1 — 9.6° (95% confidence). At p resent, there 
is no reported radio detection of this source I Kaspi et al.l 
Il996l) . however, were it one day detected, radio polari- 
metric observations could in principle constrain a. 



8 Conclusions 

The five very young pulsars with values of n measured 
via phase-coherent timing (Table [2]) show a wide range 
of spin properties and behaviors. The glitch behaviour 
exhibited by these pulsars is widely varied, ranging from 
PSR B1509-58, which has not ghtched in 21.3yr of con- 
tinuous timing observations, to the Crab pulsar, which 
experiences a glitch on average, every ~ 2yr. The mea- 
sured values of n for these five pulsars span the rel- 
atively wide range between 2.140(9) < n < 2.91(5). 
With the exception of the value of n = 2.91 ± 0.05 for 
PSR Jl 119— 6127, which is nearly compatible with n = 3, 
the measured values of n are significantly less than 3. The 
physical cause of the spin-down of pulsars remains one 
of the outstanding problems in pulsar astronomy. 
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